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Abstract 

We demonstrate control of the electron number down to the last electron in tunable few-electron 
quantum dots defined in catalytically grown InAs nanowires. Using low temperature transport 
spectroscopy in the Coulomb blockade regime we propose a simple method to directly determine 
the magnitude of the spin-orbit interaction in a two-electron artificial atom with strong spin-orbit 
coupling. Due to a large effective g- factor \g*\ = 8 ± 1 the transition from singlet S to triplet 
T + groundstate with increasing magnetic field is dominated by the Zeeman energy rather than by 
orbital effects. We find that the spin-orbit coupling mixes the T + and S states and thus induces an 
avoided crossing with magnitude A$o = 0.25 ±0.05 meV. This allows us to calculate the spin-orbit 
length Xso ~ 127 nm in such systems using a simple model. 

PACS numbers: Valid PACS appear here 
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Semiconducting nanowires are presently subject to an intense research effort due to their 
potential as nano-scale building blocks for components such as pn-iunctionSj field effect 
transistors, 3, 4| logical elements |s] and single electron circuits, js, 7^8^ 
or quantum electronic applications, e.g. qubits and quantum gates 
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10| For spintronic 



12H accurate control 

over the electron number down to the last one aswell as precise control over the spin state [3] 
and the coupling between spins in neighboring quantum dots[l^| is required. Furthermore, 
the spin-orbit interaction is of crucial importance since coupling of the spin to orbital degrees 
of freedom and thus charge determines spin relaxation and decoherence rates, which in turn 
determines decay rates for spin qubits. It is thus of fundamental interest to experimentally 
determine the spin-orbit length A50 as accurately as possible. The spin-orbit interaction 
is on the other hand also a pathway to directly manipulate the spin state of a quantum 
dot electrically rather than with conventional electron spin resonance techniques. In this 
context it was recently suggested to employ the spin-orbit coupling for all electrical spin- 



manipulation 
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17l | an approach which has the potential to be considerably faster 



than other schemes. Here, we discuss few-electron quantum dot formation in homogeneous 
InAs nanowires where the confined region is achieved by electrical depletion through local 
gate electrodes. This technique allows for highly tunable quantum dot systems with strong 
spin-orbit coupling where size, coupling and electron number starting from zero are easily 
controlled by external voltages and the spin-orbit length \$o can be determined from a 
simple measurement of the anticrossing between singlet and triplet states. 

InAs nanowires were grown using chemical beam epitaxy with 40 nm Au catalyst particles, 
yielding untapered nanowires with a diameter of ~ 50 nm. For details on growth, see Ref. Ill 
Electronic transport experiments show these wires to be n-type and have a mean free path 



lOOnm. 



19] For the fabrication of gate induced quantum dots the wires are deposited 



onto a grid of gold electrodes, which is covered by a thin SiN dielectric layer. Using this 
technique we have previously demonstrated formation of gate-defined single and double 
many-electron quantum dots. (20]] Here, due to refined sample fabrication involving thinner 
nanowires, more densely spaced electrodes (periodicity 60 nm) and a thinner SiN layer (18 
nm), the control over the potential landscape in the nanowire has been greatly improved, 
allowing tunability down to the last electron. 

The wire and five gate electrodes are contacted individually with Ni/Au contacts [see 
Fig. H](a) and inset to Fig. [11(b)] . The bias is applied symmetrically across the wire and 
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FIG. 1: (a) Sample schematic showing the Au grid covered in SiN (dark grey) and the nanowire 
(red) on top. The source and drain contacts and the five individually contacted gates are indicated, 
(b) Differential conductance measured as a function of and plunger gate voltages. The other 
gates are set to V g \ = —2.56V, Vg4 = —1.65V and V g § = 2.5V, respectively. Electron numbers are 
indicated in yellow and the inset shows an SEM image of the wire on the electrode grid, with a 1 
jj,m scalebar. 

measurements were performed in a dilution refrigerator at an electron temperature T e = 80 
mK. We apply negative voltages to gates 1 and 4 to locally deplete the electron density in the 
wire and form potential barriers, inducing a quantum dot above gates 2 and 3. Figure DJb) 
shows Coulomb blockade diamonds measured with gates 2 and 3 used as plunger gates. In 
the following we refer to the plunger gate voltage as Vq = V g3 and emphasize here that Vg 2 is 
always tuned simultaneously as in Fig. [TJb). Inside each diamond the electron number N on 
the dot is fixed, while at each diamond apex close to zero bias the dot ground state energies 
for iV and iV + 1 electrons are degenerate and transport can occur. In Fig. [0(b) no further 
degeneracy points are observed for Vq < —0.94V and the diamond borderlines continue 



without kinks, indicating that the quantum dot is empty 2l| . The diamond border-line 



conductance peaks correspond to the quantum dot chemical potential being aligned to fie 
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FIG. 2: Evolution of the differential conductance peaks as a function of B along the two cuts (a) 
and (b) in Fig. [1(b). 

or /is, as indicated in Fig. [T](b). Additional peaks in the differential conductance represent 
transitions involving excited states. 

In addition to conductance peaks originating in dot states Fig. [H(b) also shows an abun- 
dance of weak and densely spaced peaks, which are related to the low-dimensional nature of 
the source and drain leads|2l|. The unused local gates beneath the wire are set to +2.2 V 
in order to tune the leads towards metallic behavior and avoid accidental dot formation. 
However, the relatively short mean free path and a coherence length which at T=100 mK 
extends over the full nanowire, together with the radial confinement, cause a varying den- 
sity of states in the leads which leads to these resonances. This is also consistent with the 
observation of negative differential conductance lines. 

To separate these lead effects from the quantum dot states we apply a magnetic field 
B perpendicular to the substrate plane and measure the differential conductance along two 
cuts in Fig. [T](b) detecting the ground-and excited state lines of the quantum dot with 
one (a) and two (b) electrons. For the one-electron case the peaks in Fig. EJja) on either 
side of zero bias split with increasing B, reflecting transport through the Zeeman-split first 
(t, |) and second (ft, JJ-) orbital level. In low- dimensional systems the electronic g-factor g* 



depends on system size and dimensionality. It is o 



bulk modulus, which for InAs is \g*\ ~ 14.7. 
1 and 2, we extract the effective g-factors \g\ 
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ten found to be strongly reduced from the 



231 ] From the Zeeman splitting of orbitals 
8 ± 0.4 and \g%\ = 8.9 ± 1. As expected 



these values are smaller than the bulk modulus, but consistent with previous experiments 
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performed on similar nanowires. 23] 

In Fig. [2(b) we identify the two-electron ground state peak (♦) as the transition from a 
singlet S to spin-up which involves tunneling of a J,-electron and therefore moves to higher 
energy with increasing B. For the three triplet excited states T + ,T°,T~ we find two lines 
(□/■) involving transitions T + — >| and T° — >l (□) or transitions T° — >~\ and T _ — >J, (■). 
Here the first two transitions (□) require tunneling of a j-electron and the transition line 
moves to lower energy with increasing B while in the other two cases (■) a |-electron tunnels 
and the line shows a slope similar to that of the ground state line. In the following we refer 
to the two excited state lines as //yo and corresponding to the two dot states carrying 
most of the current. [24] The splitting /i T o — /i T + allows us to extract \g%,\ = 8 ±0.5 similar to 
the one-electron values. In quantum dots defined in GaAs heterostructures this splitting is 
typically only observed when a large B is applied in-plane with the 2DEG since the (/-factor 
in these systems is much smaller and orbital effects dominate for perpendicular fields. 

A quantitative analysis of our data takes into account asymmetries in the capacitive 
coupling from source and drain to the quantum dot. The conductance peak slopes in Fig. 
[1(b) change with bias and plunger gate voltage and differ between states. This indicates 
that the relation between Vq, Vd s and the energy levels of the quantum dot is non-linear. 
We therefore allow the gate leverarm «g = Cg/Cs to depend on Vq but not on Vd s and the 
source/drain coupling asymmetry 5a = as — ar> to vary with Vd s , where ats(D) — Cs(d)/Ct; 
is the source (drain) leverarm. 25|] From the two diamond border slopes at each degeneracy 



point (Vds ~ 0) we extract «g's which can be approximated as aciVc) = 0.32 — (Vq + 
0.905V)/2.27V. For symmetric coupling of the dot to source and drain it follows 5a = 0. 
Using aciVc), we re-scale Fig. [1(b) eliminating the dependence of borderline slopes on V G . 
The remaining variations in the slopes are then attributed to a V^-dependent 5a which can 
be determined from these same slopes. For non-zero 5a, the energy shift of the quantum 
dot levels with Vd s relative to source/drain is given by (1 =F Sa(Vd s )) Vd s /2 respectively. 

The necessity of this correction scheme becomes clear when we extract the single-particle 
orbital splitting Ae = 62 — £1 from Fig. d(b) and plot the bare bias splitting for negative 
(blue dots) and positive (red dots) Vd s as a function of Vq in Fig. [3(a). The splittings clearly 
deviate at most Vqs, but after the correction is applied the splittings for negative (blue 
squares) and positive (red squares) Vd s are nearly identical and we find Ae ^6.5 meV, 
changing by less than 4 % in the measured interval. The magnetic field dependence of Ae 
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FIG. 3: Orbital splitting Ae and singlet-triplet splitting Ast as a function of Vq and B. In (a) 
and (b) Ae was derived from the corrected difference between f and [see FigQ^b) and FigG(a)] . 
In (c) Ast was calculated from the data in FigQ] cut (a) and in (d) both splittings A$T+ (O) 
and A<jyo(M) were extracted from Fig[2]Jh). For comparison the small dots in panels (a) and (c) 
indicate the bare bias splittings for negative (blue) and positive (red) Vd s . 



is extracted from Fig. [2(a) and we find a quadratically decreasing behavior [see Fig. [3( 
This is similar to previous results in quantum dots with a 2D harmonic confinement. 
For comparison, we apply a similar harmonic confinement model 2l|, |28| and extract hw. 



b)l 



27] 



6.3 meV and hu y = 40 meV from the fit indicated by the black line. This is in agreement 
with the device geometry from which we expect that the dot is elongated along the nanowire. 
We further estimate an effective size 2X X = 2\Jh/m*uj x = 46 nm and 2X y = 18 nm. 

For the two-electron case we determine the singlet-triplet splitting at B = with Ast = 
Ht — Hs [see Fig. [3(c)]. Again the two bias directions give the same result after the correction 
has been applied and Ast decreases from 3.2 meV to 2.8 meV over the measured gate voltage 
range. The opposite dependence of Ast on plunger gate voltage was previously observed 
in a lateral gate-defined quantum dot in a GaAs based two dimensional electron gas and 
attributed to shape deformation of the parabolic potential. {29), [3^] 

For the magnetic field dependence of Ast [see Fig. 131(d)] we consider two splittings 
Ast+ — f^2,T+ — l^2,s an d Ast® = ^2,t° — ^2,s- ^st° only shows a small shift with B, 
which we attribute to orbital and interaction effects. Due to the large Zeeman splitting be- 
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FIG. 4: (a) Singlet-triplet transition from cut (b) in Fig. [T]after gate voltages have been optimized. 
The transition between ^+ and ^5 exhibits a clear anti-crossing, (b) Energy levels extracted from 
the corrected peak splittings in (a) together with a fit to a simple model indicated by the red lines. 

tween the triplet states, Ast+, shows a strong linear decrease with increasing B. In addition, 
we have extracted the Coulomb interaction energy U = 13.5 meV between two electrons on 
the first orbital using the corresponding conductance peak splitting in Fig. Hfb). We found 
U to be independent of Vq and B within the investigated ranges. 

Figure Ufa) again gives the magnetic field dependence of \it and \is similar to that in Fig. 
E^b) but over an extended magnetic field range. At B pa 3.8 T the quantum dot undergoes 



a transition where the two-e 
usually observed crossing 



ectro n g roundstate changes from S to T + . In contrast to the 
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321 ] . the S — T transition in our data clearly shows an 



avoided crossing due to the spin-orbit interaction. We then extract the separation of these 
three chemical potentials and by assuming that the singlet energy Es is independent of B 
we can reconstruct the spectrum of the two-electron artificial atom as shown in Fig.[4](b). 
The magnitude of the anti-crossing is A$o = 0.23 meV and serves as a direct measure of 
the spin-orbit interaction strength in such quantum dots. 



The InAs nanowires are known to have wurtzite crystal structure 23J] , despite the fact that 
bulk wurtzite InAs does not exist. The wurtzite crystal has a reduced symmetry compared 
to zinc-blende structures and because of s-p z mixing, additional spin-orbit terms linear in 



momentum occur for the conduction band electrons in bulk materials. In nanowires, like 
in most nanosystems, an important role is played by the boundary and additional spin- 
orbit terms arise in asymmetric structures, such as ours. Generally the spin couples to the 
longitudinal motion of the electron with a spin-orbit interaction that in first order of k ■ c 
assumes the following general form 

H so = (k-c)(j-a), (1) 

where c is a unit vector along the wurtzite hexagonal axis (c axis) and 7 = (7^, j y , j z ) is a 
vector of coupling constants. In a Hund-Mulliken approximation for the elongated nanowire 
quantum dot this general form of the spin-orbit interaction leads to efficient mixing of the 
S and T + states with the magnitude of the anti-crossing given by 

E z rnm| r— — E z r n ,„, 

Aso(B) = vi~v 7? + ''» = 7iw P) 

where Ez = Q*HbB is the Zeeman energy, r i2 = (i/)t\(%i — X2)\ips) is an effective distance 
between the two electrons, and ray is the electron effective mass in the nanowire in a magnetic 
field. The red lines in Fig. H^b) are a fit to the standard expression for level repulsion 



Bl , != ^ ± ^^L +A|o , (3) 

where the bare singlet energy was set to be constant Eg = and E? + = Et — g*[i^B with 
Et = —JbB. We assume that r i2 does not depend on B. This means that Aso increases 
linearly with magnetic field and yields an asymmetric splitting when comparing the left 
with the right hand side of the anticrossing. The linear dependence of Et on magnetic field 
is an experimental finding consistent with all acquired datasets and characterised by the 
phenomenological parameter 7^ = 0.14 meV/Tesla. At B w 3.8 T where the bare levels 
cross (dashed lines) we find Ez ~ 1.8 meV and estimate nm which allows us 

to calculate the spin-orbit length A50 ~ 127 nm from Eq. (J2]). This value is comparable to 



Xso ~ 200 nm obtained for unconfined electrons in similar InAs nanowires 19]. 
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Supplementary Notes 



A. Detecting the Last Electron 

B. g-factor and singlet-triplet anti-crossing for different dot configurations 

C. B-field dependence of the bare energy levels of an elongated 
quantum dot with harmonic confinement 

D. Origin of the SO-interaction in wurtzite InAs and 
avoided crossing between S and T + 
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A. Detecting the Last Electron 

To support our statement that it is possible to empty the dot completely of electrons, we 
have measured Coulomb blockade diamonds on two different samples and in different cool- 
downs. FigJ5](a) shows the differential conductance in the few-electron regime from a second 
cool-down in the dilution refrigerator using a configuration similar to the one discussed in the 
main text. Gates 2 and 3 are again used simultaneously as plunger gates, but we extended 
the sweeps to larger bias voltages in order to show the opening of the N = region without 
any features that would indicate further electrons leaving the dot. In all the measurements 
with two plunger gates the relative lever-arm between them is C^jC^ = 0.85 ± 0.05. This 
shows that the quantum dot in the nanowire is similarly coupled and tuned by both plunger 
gates and therefore likely to be centred in the middle. Figure [5](b) was measured at 4.2 K 
for the same sample. 

Figj5](c) shows Coulomb diamonds on a second sample with a similar nanowire in a 3- 




V as (mV) V sd (mV) v ds( mV ) 

FIG. 5: (a) Differential conductance for a dot with similar settings as in the main text over a larger 
bias range (V g i = —2.274V, V g 4 = —1.630V and V g § =2.2V). (b) Differential conductance for a dot 
confined by gates 1 and 3 (V g i = —1.5V, V g s = —0.7V, V g 4 and V g § are both at -0.2V) measured 
with gate 2 as plunger gate at 4.2 K. (c) Main figure: Differential conductance of a second sample 
with a similar quantum dot. Arrows indicate alignment of the chemical potential of the dot with 
that of the drain. Inset: absolute value of the current trough the same device. 
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gate configuration. Here, the last electron leaves the quantum dot at V g = —0.3 V. This is 
supported by the fact that no further borderlines (see black arrows) appear at positive bias 
after the chemical potential of the drain contact has dropped below the chemical potential of 
the last electron on the dot (lowest arrow). In other samples [see main text or Figf5](a) and 
(b) ] similar but less distinct lines (both differential and negative differential conductance 
lines) do occur even below the N=l diamond. This is something that we observe in many of 
our nanowire based quantum dots. Specifically, also in the InP/InAs heterostructure based 
QD's where the electron number can be determined without any doubt [23]. The origin of 
these peaks lies in the fact that our leads to the QD's are nearly one- dimensional (in contrast 
to 2DEG based dots). A relatively short mean free path and a coherence length at T = 100 
mK which extends over the full nanowire, lead to the observation of resonances that have 
their origin in the leads and not the quantum dot. Such an explanation is also consistent 
with the observation of negative differential conductance lines. We use the magnetic field 
dependence of the differential conductance peaks to determine which peaks belong to the 
dot and which ones originate in the leads. The dot states typically exhibit a " slow" variation 
as a function of magnetic field while the leads states vary quickly in energy with B and tend 
to be suppressed for large fields [ see e.g. Fig. 2 (b), (c) and Fig. 3 (b), (c) ]. In addition, 
by tuning the gates which lie below the contacts (left of gate 1 and right of gate 4) we were 
able to tune the lead states while leaving the dot related peaks intact, therefore excluding 
the possibility that these additional resonance lines originate from low lying excited states 
in the dot. 

We have fabricated and performed measurements at low temperatures on six devices with 
a gate period of 60 nm. Two of these had five or more working gates without any leakage or 
connection problems and in both we were able to access the last electron in a gate defined 
quantum dot. 
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B. g-factor and singlet-triplet anti-crossing for different dot configurations 

The current device layout permits the definition of very flexible quantum dots, where 
couplings and the quantum dot position and geometry can be varied. In the main text, 
we consider a fairly extended dot, tuned by double plunger gates (4-gate dot). Due to the 
tunability of our device design, we can also define smaller dots only with a single plunger 
gate (3-gate dot). As the quantum dot is tuned between these two configurations, with either 
one or two plunger gates, the main features of the differential conductance as a function of 
bias and plunger gate voltage are preserved. Most notably the avoided crossing between 
the singlet and triplet two-electron states is found for both gate settings. FigsEl and [7] 
show similar measurements as Figs. 1 and 2 in the main text, but for two additional gate 
configurations. 

Fig. [6] shows Coulomb diamonds and magnetic field dependence of the conductance peaks 
for a three gate dot with gate 2 acting as the plunger gate. The Zeeman splitting and first 
orbital excited state as well as the triplet states are clearly identifiable. From the Zeeman 
splitting of the first and second orbital states in FigJH^b), we extract effective g-factors 




v d5 {mV) V dE (mV» 



FIG. 6: (a) Differential conductance measured for a quantum dot confined by gates 1 (V^i = 
—2.56V) and 3 (V g s = —1.72V). Gates 4 and 5 are both at 2.5 V and gate 2 is used as the plunger 
gate, (b) and (c) Evolution in magnetic field of the conductance peaks along cuts (b) and (c) in 
panel (a). 
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\g*\ = 8.6 ± 0.3 and {g^l — 7.6 ± 0.5, which are again in the same range as for the dot in 
the main text. The discrepancy in the magnitude of the Zeeman splitting (clearly visible 
in FigJ^b)) is reduced to a negligible difference when the correction scheme described in 
the main text is applied. In Fig. [6](c) we can again identify the transitions from the single- 
electron to a singlet or triplet state. At B= 2.5T the peaks anti-cross. Figure [7] shows the 
corresponding data for a quantum dot defined at intermediate gate settings. Compared to 
Fig. El V 9 3 is tuned more positive and V g 4 more negative, resulting in a configuration where 
the right-hand barrier is shared between gates 3 and 4. From the splitting of the states in 
Fig. Ulb) we can again extract effective g-factors which are found to be \gl\ = 7.7 ± 0.5 and 
1^1 = 7.6 ± 0.5. 

The cut for the two- electron case in FigJTKc), shows a charge rearrangement at B = 1.5 T. 
Extracting the singlet and triplet energies (see inset) shows, however, that there has been no 
change in the relative positions of the conductance peaks. Again the anti-crossing is clearly 
present in this case. 



, » 4 8 

(a) ci=^a^ta G(io 3 xe 2 /h> {b} 




V d5 (mV) V ds (mV) 

FIG. 7: (a) Differential conductance measured for a quantum dot confined by gates 1 (V g i = 
-2.56V) and 3/4 (V g3 = -1.17V and V g4 = -0.5V). V g5 = 2.5V and gate 2 is used as the sole 
plunger gate, (b) and (c) Evolution in magnetic field of the conductance peaks along cuts (b) and 
(c) in panel (a). The inset in (c) shows the singlet and triplet energies and the anticrossing at 
B m 3.5 T. 
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C. B-field dependence of the bare energy levels of an elongated quantum dot with 
harmonic confinement. 



The B-field dependence of the bare energy levels (without spin-orbit interaction) of a 
quantum dot with harmonic confinement in all three directions can be calculated by diag- 
onalizing a quadratic form 43J. For the weakest confinement direction along the axis of 
the wire (x-axis) and a transverse magnetic field B\\z, the lowest excitation energy of a 
single-electron quantum dot is given by 



A£ =2 



UJ V + Ll) x Y + LOi 



(4) 



where hu x and hio y are, respectively, the confinement energies for motion along the wire 



and in the direction transverse to both wire and magnetic field; ui c = eB / c^Jmpn^ is the 
cyclotron frequency. For small u> c ,u> x <C u y , one obtains a quadratic B-field dependence of 
the excitation energy 



Ae = huj 7 



2lu 2 

V 



(5) 



The values hu x = 6.3 meV and hw y = 40 meV given in the main text were determined using 
a fit to the measured data points Ae(B) in Fig. 3(b) of the main text. 
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D. Origin of the SO-interaction in wurtzite InAs and avoided crossing between S 
and T+ 



The spin-orbit interaction in Eq. (1) of the main text is peculiar to nanowires in which 
the transverse quantization is much stronger than the longitudinal one. In the absence of 
magnetic fields, the lowest electron mode in the nanowire is doubly degenerate, due to the 
time-reversal symmetry, and thus the dispersion relation obeys 

= e i ( 6 ) 

where e s (k x ) is the energy of the electron with longitudinal momentum k x and spin s quan- 
tized along the direction of the internal field. The spin-splitting of the mode is absent at 
k x = and can be expanded in powers of k x 

e,fe)-e^) = ^ + em p) 

In Eq. (1) of the main text, we retain, for simplicity, only the linear in k x terms. Inclusion 
of an arbitrary odd function of k x does not change our results substantially, while it can be 
viewed as redefining the 7's in Eq. (2) of the main text. 



The linear in k terms in Eq. (1) of t 



interaction of the bulk wurtzite 



C[k x c] • cr, where 

c is a unit vector along the c-axis of the wurtzite (uniaxial) crystal. In the case of wurtzite 
InAs nanowires, c is also the unit vector along the wire. Unlike in zinc-blende alloys, the 



re main text are different from the Rashba spin-orbit 



(i) 



Rashba terms are present in wurtzite in the bu 



33 



k and values of C ~ (0.1 - 10) x 10" 10 eV cm 



36 



371 ] . For wurtzite InAs, the value of C is 



have been reported for AnB V i compounds 
not known. 

Despite the fact that the Rashba terms, if present, might couple the spin efficiently to the 
transverse electron motion, these terms do not mix singlet and triplet states in nanowires 



with soft confinement along c. Notably, after averaging over the transverse electron motion 
the Rashba terms vanish. We have further verified that H^q gives no contribution to A$o 
up to the third order of perturbation theory in HgQ, provided the quantum dot electron has 
well separated scales of a tight transverse and extended longitudinal motion. In the opposite 
limiting case (not realized in our experiment), when the gate-defined confinement along the 
wire is stronger than the transverse quantization, we expect Aso to be a measure of H^q. 
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Similar to zinc-blende compounds, third and higher (odd) powers of k also contribute 

(2) 

to the spin-splitting in the conduction band of wurtzite. We denote these terms by Hg/^ 



h(k)-cr, where h(k) obeys h(kc) = 0, Vfc, due to the uniaxial symmetry of the wurtzite 35]. 
At the third order of k, h(k) reads [38] h(k) = (Aj(fe • c) 2 + X t (k x c) 2 ) [k x c], where A; and 
At are coupling constants. At higher orders in k, we expect h(k) to give a finite contribution 
to Eq. (1) of the main text, after averaging over the transverse electron motion. 

Apart from the crystalline inversion asymmetry, an important source of spin-orbit inter- 
action in nanowires is the structure inversion asymmetry at the radial edge of the sample. 
This mechanism is analogous to the Bychkov-Rashba spin-orbit interaction in heterostruc- 



tures [42J and in quantum wires it gives a finite contribution to Eq. (1), provided the electron 
is displaced from the center of the wire to one of its edges. Due to the particular way of 
gating our quantum dots (with gates beneath the wire) we expect a strong contribution from 
this mechanism in our specific case. Additional possible sources of spin-orbit interaction in 
nanowires include built-in radial strain, charged surface states, deviations of the InAs lattice 
from wurtzite (e.g. occurrences of the zinc-blende stacking sequence along c), etc. 

Next we show that the spin-orbit interaction in Eq. (1) of the main text mixes the singlet 
(S) and triplet (T + ) states and derive Eq. (2) of the main text. Implying the Hund-Mulliken 
approximation, we write the singlet and triplet orbital wave functions as follows 



ips(ri, r 2 ) = . 1 j9 (^+(ri)^+(r 2 ) - </>V>-(ri)^-(r 2 )) , (8) 
V 1 + 



iPt(vi, r 2 ) = -j= (^+(ri)V>-(r 2 ) - ^_(ri)^ + (r 2 )) , (9) 



where <j) is the so-called interaction parameter 



39 



40 



4l| and if> ± (r) = exp(iM)f±(x)g(y, z), 



with f±(x) = ±/-t-(— x) being the trial wave functions that minimize the energy of the 
two interacting electrons in the quantum dot. Here, we choose x parallel to the nanowire 
and z along the magnetic field. The orbital effect of the magnetic field can be accounted 
for by the operator M, which assumes M = ^ 2 J^. — h 0(B 2 ) in the Landau gauge 
A(r) = (—yB, 0, 0). M mixes the wave function components g(y, z) and f±(x) to an amount 
controlled by the transverse oscillator frequency uo y . However, in our further derivation the 
orbital magnetic field is not essential and we retain it only for the sake of completeness. 
Considering further a generic Zeeman field, 

H z = l -E z n ■ a, (10) 



where Ez is the Zeeman splitting and n is the spin quantization axis, we transform away 
the spin-orbit interaction obtaining 

H S o + Hz^Hz + Hl°, (11) 

where H§° = E z {myh 2 )(r ■ c) [ n x 7] ■ cr gives the combined effect the spin-orbit and 
Zeeman interactions. The longitudinal effective mass is slightly larger due to the orbital 
magnetic field effect, 

1 1 r . 2 

(12) 



1 1 

rhn ml 



1- ^ 



2wJ + w? 



However, using ftu y = 40 meV and for B < 5 Tesla this correction is negligible and we 
assume rn| = mjj = = m£ uifc = 0.023 m e for all numerical estimates. 

Taking matrix elements of if§° for two electrons and between the singlet and triplet 
states, we obtain the splitting 2A so , with 

A so = |<* 5 |tff°|* T+ >| = Eznz-jJL^-ij-ny. (13) 

Next we evaluate the distance r i2 = (i[>t\(xi — X2)\i>s) i n t wo limiting cases: (i) strong 
Coulomb interaction and (ii) weak Coulomb interaction. For our quantum dots, case (ii) 
turns out to be relevant. 

For strong Coulomb repulsion, we construct f±(x) from two Gaussian wave packets cen- 
tered around x = ±d/2 and with overlap S. Then we obtain 

1 + d 
v /2(l + 02) v ^^ 1 } 

A starting approximation for d is found by minimizing the potential energy of the two 
electrons, W = e 2 / ' nd + m*^oj 2 x d 2 /A. The result reads d = \ x (2\ x /a B ) 1/3 « 25 nm, where 
a* B = h 2 K,/m^e 2 is the Bohr radius, with k 15 being the dielectric constant of bulk InAs, 
and X x = ^Jh/m^u x = 23 nm the effective longitudinal dot size. Similarly, the width of 
the Gaussian wave packet Ao can be estimated by minimizing the sum of the confining and 
exchange energies, S = h 2 /2mt\l + S 2 e 2 /nd ) where S = exp(— d 2 /A\l). The result for 
Ao is then given by Aq = d [2 ln^/a^)] -1 ^ 2 , where d* B = h 2 K/mte 2 is the effective Bohr 



radius. Similarly, we find S = ^a* B /d. In the strong Coulomb interaction regime, defined 
by \ x /a* B ^> 1 and d/d* B ^> 1, one can further evaluate with the help of the expression 
= y/l + (At H /U H ) 2 - At H /U H and by letting t H ~ Sh 2 /m\\l and U H ~ e 2 / k\ . 
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In the absence of Coulomb interaction, we obtain 



h 



r\2 = 



1 



cut. 



1/4 



(15) 



We note that this expression gives a lower bound for r i2 in the presence of Coulomb inter- 
action. 
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